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RADIATION  PATTERNS  OF  A  DUAL  REFLECTOR  SYSTEM 
BASED  ON  THE  GEOMETRICAL  OPTICS  APPROXIMATION 

1.  INTRODUCTION 

In  a  previous  report  {1],  an  electronically  steerable  dual  reflector  sys¬ 
tem  fed  by  a  phased  array  was  discussed.  The  configuration  of  such  a  systems 
in  shown  in  Figure  1*.  The  cross  section  of  the  main  reflector  of  this 
system  is  a  parabola  while  the  subreflector  is  an  ellipse.  The  novelty  of 
this  system  is  that  the  two  foci  of  the  ellipse  are  located  respectively  at 
the  center  of  the  main  reflector  and  the  center  of  the  feed  array.  There¬ 
fore,  rays  from  all  directions  reflected  at  the  center  of  the  main  reflector 
will  converge  at  the  center  of  the  feed  array.  If  the  scan  angle  is  small 
and  the  sizes  and  locations  of  both  reflectors  are  properly  chosen,  rays 
within  the  required  scan  range  will  be  distributed  on  the  array  face.  The 
required  feed  array  and  subreflector  sizes  are  minimized  and  the  spillover 
loss  can  be  very  small. 

In  reference  [1],  a  method  for  designing  such  a  system  is  discussed,  which 
is  based  on  the  geometrical -optics  method.  The  same  approach  will  be  used  in 
this  report  to  compute  a  two-dimensional  radiation  pattern  of  such  a  system. 
It  is  known  that  the  geometrical  optics  approach  is  only  an  approximation, 
particularly  when  interference  effects  are  sought.  To  accurately  compute  the 
radiation  patterns  of  such  a  system,  a  current  distribution  integral  method 
must  be  used.  However,  such  an  approach  is  tedious  and  requires  a  major  ef¬ 
fort.  At  the  present  time,  we  are  only  interested  in  a  preliminary  study  to 
determine  the  feasibility  of  such  a  system.  Therefore  the  geometrical-optics 
method  is  used.  If  the  results  are  promising,  effort  will  be  applied  to  the 
current  distribution-integral  method. 

In  this  approach,  a  receiving  case  is  assumed.  All  Incoming  waves  are 
assumed  to  be  planar  with  various  incidence  angles  with  respect  to  the  main 
reflector.  At  each  Incidence  angle  a  large  number  of  rays  is  assumed  and  they 
are  uniformly  distributed  on  the  main  reflector.  The  path  length  of  each  ray 
from  the  incident  wave  front  to  the  feed  array  is  then  computed,  and  the  number 
of  rays  Intersecting  each  element  in  the  feed  array  is  determined.  Thus,  the 
amplitudes  and  phases  of  these  incoming  waves  on  the  feed  array  can  be  deter¬ 
mined,  and  the  radiation  pattern  can  then  be  computed. 

In  this  report,  we  wish  to  determine  whether  it  is  feasible  to  steer  the 
beam  off  the  broadside  and  how  this  will  affect  the  sldelobe  level.  Secondly, 
wa  would  like  to  determine  whether  a  conventional  method  can  be  used  for  pat¬ 
tern  synthesis.  Third,  we  wish  to  know  how  closely  the  feed  elements  in  the 
feed  array  must  be  spaced  and  how  they  can  be  phased  to  achieve  beam  steer¬ 
ing. 

2.  PATTERN  COMPUTATION  BY  USE  OF  GEOMETRICAL-OPTICS  METHOD 

Figure  1  shows  a  double  reflector  system  where  the  surface  R  is  the  main 
reflector  and  S  is  the  subreflector.  The  feed  array  is  labeled  A.  Assume  that 
a  plane  wave  propagates  towards  the  main  reflector  R  with  an  Incidence  angle 
8.  To  compute  the  radiation  intensity  at  each  incidence  angle  with  respect  to 
the  main  reflector,  it  is  assumed  that  a  large  number  of  parallel  rays  are 
incident  on  the  main  reflector  and  that  they  are  uniformly  distributed  in  a 


plane  perpendicular  to  their  direction.  These  rays  are  reflected  In  accord¬ 
ance  with  Snell's  law  by  both  the  main  reflector  and  the  subreflector,  and 
they  finally  intersect  the  feed  array.  To  compute  the  radiation  pattern  the 
phase  and  amplitude  of  the  incident  wave  at  the  feed  array  must  be  known. 
At  each  incidence  angle,  the  ray  phase  is  determined  by  the  distance  traveled 
by  the  ray  from  the  plane  wave  front  to  the  main  reflector  and  the  doubly 
reflected  path  from  the  main  reflector  to  the  array  surface.  The  ray  path 
length  from  the  wave  front  to  a  point  on  the  main  reflector  (£q)  is  computed 
by  using  equation  (A3)  given  in  Appendix  A.  The  computations  of  the  reflected 
ray  path  lengths  from  the  main  reflector  to  the  subreflector  (Aj)  and  from 
subreflector  to  the  feed  array  ( I2 )  are  given  in  Appendix  B.  The  phase  of  each 
ray  at  the  array  face  is  then  modulo  2n  of  (*  /X)  (Aq+Aj+A?)  where  X 
is  the  wavelength  of  the  incident  wave.  Since  both  main  reflector  and 
subreflector  are  curved  surfaces,  the  rays  that  are  uniformly  incident  on  the 
main  reflector  are  not  necessarily  uniformly  Incident  on  the  feed  array  face. 
The  number  of  rays,  or  flux  incident  per  unit  area  at  a  point,  represents  the 
intensity  of  the  incident  field  at  that  point.  If  a  large  number  of  rays  (or 
flux)  is  assumed  incident  on  the  main  reflector,  and  the  feed  array  is  divided 
into  n  sections  (or  n  elements)  where  n  is  a  much  smaller  number  than  that  of 
the  total  number  of  rays,  then  the  number  of  rays  incident  on  each  of  the  n 
sections  represents  the  intensity  of  the  incident  field  on  that  element.  For 
pattern  synthesis,  a  weighting  function  is  usually  applied  to  each  element  to 
achieve  the  desired  low  sidelobe  level.  For  the  present  case,  this  weighting 
function  must  take  into  account  this  non-uniformity  property.  A  critical 
design  problem  relates  to  the  variation  of  incident  field  intensity  on  each 
element  as  a  function  of  scan  angle.  If  it  changes  too  much  at  different 
scan  angles,  it  may  be  impossible  to  achieve  a  desired  pattern  with  a  constant 
array  weighting  function.  However,  the  examples  which  we  have  studied  indi¬ 
cate  that  the  field  intensity  on  the  array  face  is  uniformly  distributed  and 
its  pattern  resembles  that  of  a  conventional  linear  array.  Of  course,  in 
arriving  at  this  conclusion,  it  is  assumed  that  the  array  element  is  not  a 
point  source,  has  omnidirectional  element  pattern,  and  only  receives  energy 
from  the  direction  of  the  wave  front,  or  in  the  ray  tracing  region. 

To  steer  the  beam  to  an  angle  8q,  the  required  phase  setting  at  each  ele¬ 
ment  can  be  achieved  by  averaging  the  path  lengths  of  all  rays  which  intersect 
that  element  and  which  are  incident  at  an  angle  0q  to  the  main  reflector.  The 
required  element  phase  setting  is  just  the  conjugate  of  this  average  phase. 
One  may  see  that  if  the  path  phase  of  each  individual  ray  within  each  element 
deviates  too  much  from  the  average  phase  value,  then  it  may  be  difficult  to 
form  a  beam  in  the  9q  direction.  This  can  happen  when  the  array  element  spac¬ 
ing  is  too  large. 

3.  COMPUTED  RESULTS 

As  an  example,  the  radiation  pattern  of  a  double  reflector  system  is 
computed  and  some  of  its  results  are  presented  here.  The  system  computed 
here  is  designed  by  the  method  discussed  in  reference  [1].  The  dimensions  of 
this  system  are  as  follows: 

Main  Reflector  50X 


Sub  Reflector 


20X 


Feed  Array  16X 

Focal  Length  (Main  Reflector)  84X 

Some  computed  patterns  for  scan  angles  at  -2  degrees,  0  degree  and  8 
degrees  are  shown  respectively  in  Figures  (2a),  (2b)  and  (2c).  The  feed 
array  is  uniformly  excited  and  has  an  element  spacing  of  X/2.  The  patterns 
are  plotted  in  a  range  from  -20  degrees  to  20  degrees  which  is  the  required 
field  range  of  the  system  contemplated. 

Similar  patterns  are  plotted  for  the  same  system  in  Figures  (3a),  (3b), 
and  (3c)  at  similar  scan  angles.  However,  in  this  case,  the  feed  array  is 
excited  by  a  30  dB  Chebyshev  weighting  function. 

Figures  (4a)  and  (4b)  show  the  radiation  pattern  of  the  same  system 
except  that  the  element  spacing  is  .9X  Instead  of  .5X. 

From  these  figures  the  following  conclusions  are  drawn: 

a.  When  the  feed  array  is  excited  uniformly,  its  radiation  pattern  is 
similar  to  that  of  a  conventional  uniformly  weighted  linear  array.  No 
appreciable  degradation  is  evident  when  the  beam  is  steered  over  a  range  from 
-2  degrees  to  8  degrees. 

b.  A  conventional  weighting  function  such  as  the  Chebyshev  designed  for 
30  dB  sldelobe  yields  similar  results  to  those  obtained  with  a  conventional 
linear  array.  However,  when  attempts  were  made  to  further  reduce  sidelobes, 
they  do  not  seem  to  work.  This  will  be  discussed  in  detail  in  the  next 
section. 

c.  When  element  spacing  Increases  beyond  a  certain  limit,  a  kind  of 
grating  lobe  appears.  However,  the  exact  point  of  limitation  is  not  clear  at 
this  point.  It  varies  from  system  to  system.  Further  study  in  this  area  is 
required. 

4.  SIDELOBE  DESIGN 

Since  the  achievement  of  low  sidelobes  is  very  important,  attempts  were 
made  to  achieve  lower  sldelobe  levels  than  in  the  examples  shown  in  section 
3.  The  first  attempt  is  to  use  a  conventional  array  approach  in  which  40  dB, 
50  dB  and  60  dB  Chebyshev  array  weights  are  used.  Patterns  for  such  illumina¬ 
tion  weights  are  shown  in  Figures  (5a),  (5b)  and  (5c),  respectively.  The 
configuration  of  the  system  of  these  plots  is  the  same  as  that  of  Figures  2 
and  3.  One  may  see  from  these  figures  that  the  attainable  sldelobe  level 
seems  to  be  limited.  Another  approach  to  tapering  the  array  illumination  was 
also  tried,  as  shown  in  Figure  6.  This  array  illumination  function  has  a 
cosine  squared  distribution  with  edge  element  20  dB  down  from  the  center 
element.  Again,  the  achievable  sldelobe  level  is  quite  limited.  Patterns  at 
other  scan  angles  for  both  illumination  cases  are  similar;  hence  they  are  not 
shown. 

In  order  to  investigate  why  a  conventional  pattern  synthesis  method  does 
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not  work  in  this  system,  the  phase  and  amplitude  distribution  across  the 
array  surface  for  plane  waves  incident  at  0  and  -2  degrees  to  the  main 
reflector  are  plotted  in  Figures  (7a)  and  (7b).  Figure  (7a)  shows  the 
amplitude  distribution  (or  intensity  of  incident  waves).  Although  it  varies, 
it  is  reasonably  uniform.  Compensation  of  this  variation  with  array  illumina¬ 
tion  function  does  not  seem  to  yield  the  desired  result.  Amplitude  distribu¬ 
tions  of  plane  waves  incident  at  other  angles  show  similar  properties. 

Figure' (7b)"  shows  -t-he- -phase— di-atribntion  for__the  cases  of  0  and  -2  de¬ 
grees.  (Actual  Phase,  not  mod.  2m )  This  phase  distribution  at  each  -  - 
element  is  actually  the  average  of  the  phases  of  all  the  incident  waves 
which  intersect  with  that  array  element.  One  may  see  that  the  variation  of 
the  phase  across  the  array  face  has  a  quadratic  distribution.  A  conventional 
linear  array  has  a  linearly  progressive  phase  distribution  and  its  pattern 
synthesis  is  based  on  this  property.  It  is  evident  that,  for  an  array  of 
this  type,  a  conventional  linear  array  pattern  synthesis  method  cannot  pro¬ 
duce  the  desired  result.  Further  work  along  this  line  is  required. 

5.  CONCLUSIONS 

The  results  presented  in  this  report  are  encouraging.  The  radiation 
pattern  behave  reasonably  well.  There  seems  to  be  no  difficulty  in  beam 
steering  and  it  also  responds  well  to  the  conventional  pattern  control 
weighting  function  to  a  certain  degree.  However,  there  are  certain  problems 
which  require  further  investigation,  such  as  array  element  spacing  and  pattern 
synthesis,  it  has  been  pointed  out  that  since  these  results  are  based  on  the 
geometrical-optics  method,  it  is  only  an  approximation.  Because  of  these 
encouraging  results,  we  recommend  that  further  effort  be  extended.  We 
recommend  that  this  problem  be  rigorously  formulated  by  use  of  the  current- 
distribution  integration  method,  and  we  also  recommend  that  the  three 
dimensional  case  be  analyzed. 
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Fig.  2(b)  —  Radiation  pattern,  steering  angle  =  -2°,  element  spacing  =  X/2. 


Radiation  pattern,  steering  angle  =  8°,  element  spacing  =  A/2 


Radiation  pattern,  steering  angle  =  0°,  element  spacing  =  X/2,  30  dB,  Chebyshev  weighting. 


Fig.  3(b)  —  Radiation  pattern,  steering  angle  =  -2°,  element  spacing  =  1/2X,  30  dB  Chebyshev  weighting. 


Fig.  3(c)  —  Radiation  pattern,  steering  angle  =  8  ,  element  spacing  =  X/2,  30  dB  Chebyshev  weighting. 


Radiation  pattern,  steering  angle  =  0U,  element  spacing  =  ,9X,  30  dB  Chebyshev  weighting. 


Fig.  4(b)  —  Radiation  pattern,  steering  angle  =  -8°,  element  spacing  =  ,9X,  30  dB  Chebyshev  weighting. 


Fig.  6(b)  —  Radiation  pattern,  steering  angle  =  0°,  element  spacing  =  X/2,  30  dB  Chebyshev  weighting. 


Radiation  pattern,  steering  angle  =  0U,  element  spacing  =  A/2,  60  dB  Chebyshev  weighting. 


Radiation  pattern,  steering  angle  =  0°,  element  spacing  =  A/2,  cosine  square  tapering  weight. 
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Pig.  7(a)  —  Array  amplitude  distribution  of  plane  wave  incident  at  0 
and  -2°  with  respect  to  the  main  reflector. 


APPENDIX  A 


Figure  AI  shows  a  wave  front  of  a  plane  propagate  towards  the  main 
reflector  R  with  an  incident  angle  0  .  From  this  figure, 

rp  -  (Xn-Xf)  (Al) 

pq  -  RqP  tanQ  (A2) 

-  (Ym-Yf)tan0 

Uol  -  (rp+pq)cos9 

-  [ Xm-Xf )+( Ym-Yf )tan0 ]cos0  (A3) 

where  (X^Yg,)  is  a  reference  point  and  (Xf,Yf)  is  a  point  on  the  main 
reflector  It  can  be  shown  that  this  equation  also  applies  to  the  case 
when  9  is  negative. 


\ 
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Fig.  A-l  —  Path  length  from  wave  front  to  main  reflector 


APPENDIX  B 

Let  a  parabola  R  be  described  by  an  equation; 

y2  -  4fa(x+fa)  (Bl) 

where  fa  Is  focal  length  of  the  parabola. 

The  gradient  of  such  a  curve  Is  then 

7  G  -  -2fm  1  +  yj  (B2) 

where  1  and  j  are  respectively  unit  vector  along  x  axis  and  y  axis.  The  unit 
normal  vector  A  is  then 

n  -  2fi  -  yi 

n — r  •  (B3) 

Jy2+4fa2 

An  incident  wave  having  Incident  angle  0  with  the  x-axis,  can  be 
represented  (See  Figure  Bl) . 

Zq  -  cos0  i  +  sin0  j  . 

The  reflected  wave  can  be  represented 

£I  "  ~l0  +  •!<>> 


x. cos0+y sin 6  ^  x. sin0-y. cos0. 

■  -4f*  2.,.2  ««.  2... 2  1  J 


where  (x^yj)  is  the  point  on  the  main  reflector  where  the  ray  Zi  is  reflected. 
This  Zj  can  be  written  into  an  equation  of  straight  line  which  has  the  form 


y  -  mfx+kf 


(B5a) 


(y12-Kt12)tan0 


(B5b) 


y.-Xjtan©  (y.  )tan0 

'h,r*  f  *  ,7iw1  kf  ■  ,^8^ —  •  <B5b> 

In  order  to  find  the  length  Zj  of  the  reflected  wave  propagating  from 
the  main  reflector  to  the  subreflector,  it  is  required  to  know  the  equation 
which  describes  the  elliptical  subreflector.  This  equation  is 


where 


a^2-tayy2+(bx+2ty)x+byy+cxy  -  0 


ax  "  ^e-V2-4*2 

*y  -  <ye-yd)2*4a2 


(B7a) 

(B7b) 
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bx  -  <xe_xd>axy+8a2xe 

(B7b) 

(B7c) 

by  -  (ye-yd)axy+8a2ye 

(B7d) 

t  -  (Xe-xd)(ye-yd) 

(B7e) 

Sy  ■  (xd2+yd2)"(xe2+ye2)'4a2 

(B7f  ) 

where  (xg,ye)  and  (xd,yd)  are  respectively  the  two  focal  points  of  the  ellipse 
located  at  the  center  of  the  main  reflector  and  the  center  of  the  feed  array 
and  a  is  the  distance  between  these  two  focal  points.  Since  the  ellipse  is 
assumeed  to  pass  through  the  origin ,  cXy  *  0. 

The  intersection  point  of  the  reflected  ray  and  the  subreflector  is 

then 


The  sign  is  so  chosen  to  yield  the  larger  value  of  xs ,  where 


(B8) 


b  +2tk.+2a  mrk-+b  m£ 

2B - * - y  1 

a  +2  tnua  m. 
x  f  y  f 


(B9a) 


a  k£  +b  k 

c  -  -xi — i-L 


ic+2 1  f  mj +aymf 4 


(B9b) 


mf  and  kf  are  parameters  of  the  reflected  ray  Jlj  which  are  defined  in 
equation  (B5b),  and  bx,  t,  ay,  by  are  defined  in  equation  (B7a)  through  (B7f). 

The  y  coordinate  of  this  point  can  be  found  from  (Eq.  B5a)  when  x  is 
known.  The  length  of  is  then 
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I  1 1 1  -  txs-xi >2  +  (yg-y^2!  • 

The  normal  vector  n'  of  the  ellipse  is 
o'  1  -*  n'  o 

«'-===J-  +  7===.  j 

^n'  ^+n*  yj  n'  ^+n'  2" 

where 

n'l  *  2xgax+2tyg+bx 

n’2  *  2ayyg+2txg+by 

where  (ys,xs)  is  a  point  on  the  subreflector. 
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(B12a) 

(B12b) 
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